Cyclotron Maser Emission from Blazar Jets? 
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ABSTRACT 



We consider the production of electron cyclotron maser emission by low- 
density, highly magnetized plasmas in relativistic jets. The population inversion 
required to drive cyclotron maser instability could occur in localized, transient 
sites where hydromagnetic instabilities, shocks, and/or turbulence lead to mag- 
netic mirroring along current-carrying flux tubes. The maser is pumped by the 
conversion of kinetic and magnetic energy into j • E work as electrons are accel- 
erated into the maser by the parallel electric field that develops as a result of 
the mirror. We estimate the maximum brightness temperatures Tb that can be 
obtained in a single maser site and in an array of many masers operating simulta- 
neously, under conditions likely to apply in blazar jets. Synchrotron absorption, 
by relativistic electrons within the jet, presents the largest obstacle to the es- 
cape of the maser radiation, and may render most of it invisible. However, we 
argue that a high brightness temperature could be produced in a thin boundary 
layer outside the synchrotron photosphere, perhaps in the shear layer along the 
wall of the jet. Induced Compton scattering provides additional constraints on 
the maximum brightness temperature of a masing jet. We suggest that recent 
observations of diffractive scintillation in the blazar J1819+3845, indicating in- 
trinsic brightness temperatures greater than 10 14 K at 5 GHz, may be explained 
in terms of cyclotron maser emission. High- Tb maser emission from blazar jets 
may extend to frequencies as high as ~ 100 GHz, with the maximum possible 
Tb scaling with frequency roughly oc v^ 1 . Less massive relativistic jet sources, 
such as microquasars, are even better candidates for producing cyclotron maser 
emission, primarily in the infrared and optical bands. 
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Subject headings: acceleration of particles — masers — radiation mechanisms: nonthermal— 
galaxies: jets — radio continuum: galaxies 

1. Introduction 

There has been a longstanding debate over whether some compact extragalactic radio 
sources have intrinsic brightness temperatures that exceed the maximum value consistent 
with incoherent synchrotron radiation. The maximum brightness temperature for steady 
emission, ~ 10 12 K, is set by the " inverse Compton catastrophe" (Kellermann & Pauliny-Toth 
1969), although attainment of this limit requires very large deviations from equipartition 
between magnetic and particle energy and, therefore, unacceptably high energy requirements. 
A more plausible upper limit, ~ 10 11 K, is obtained for self-absorbed synchrotron sources 
close to equipartition; most resolved compact sources fall near or below this limit (Readhead 
1994). 

However, the observational limits on brightness temperature are not as clear-cut for 
bright, unresolved radio sources that exhibit rapid "intraday variability" (IDV) at GHz 
frequencies. If simple light-travel time arguments are used to infer an angular size from 
the variability timescale, then some of these sources must have brightness temperatures as 
high as T B ~ 10 21 K. The actual angular size of the source can be a factor ~ V larger than 
the naive estimate if the radiating material (or, rather, its emission pattern) approaches 
us with a bulk Lorentz factor T ^> 1. The solid angle is thus a factor ~ T 2 larger, while 
the intrinsic brightness temperature of the beam is boosted by a factor V relative to the 
comoving frame of the source. Thus, the apparent Tg could be boosted by a factor ~ T 3 
over its value in the comoving frame, and values of V ~ 10 3 could explain most extreme 
instances of IDV. However, such large bulk Lorentz factors also imply very low synchrotron 
radiative efficiencies, leading to implausibly large jet energy fluxes (Begelman, Rees & Sikora 
1994). 

The need for high intrinsic brightness temperatures appeared to become less acute when 
it was confirmed that IDV results mainly from refractive interstellar scintillation (Kedziora- 
Chudczer et al. 1997; Dennett-Thorpe & de Bruyn 2000, 2002). These observations do not 
require intrinsic brightness temperatures larger than ~ 10 13 K, which are easily explained by 
incoherent synchrotron radiation boosted by a modest Lorentz factor < 30. However, space 
VLBI has failed to resolve these and many other compact sources (Tingay et al. 2001), 
leaving open the possibility that the brightness temperatures are much higher. Moreover, 
recent observational evidence for diffractive scintillation from the blazar J1819+3845 (Mac- 
quart & de Bruyn 2004) yields an upper limit for the angular size of < 10 yuas, implying 
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intrinsic brightness temperatures in excess of 10 14 K. To explain such high intrinsic Tb with 
an incoherent synchrotron model would require Y > 10 3 , once again leading to implausible 
energy requirements. 

Brightness temperatures that exceed the synchrotron limit, after relativistic effects are 
taken into account, could indicate that a coherent radiation mechanism is at work. There 
have been several attempts to apply such mechanisms to the conditions in relativistic jets 
from active galactic nuclei (e.g., Baker et al. 1988; Sol, Pelletier & Asseo 1989; Krishnan 
& Wiita 1990; Lesch & Pohl 1992; Benford 1992; Benford & Tzach 2000). All of these 
models rely on the development of strong Langmuir turbulence, and implicitly assume that 
the electron plasma frequency v v exceeds the electron cyclotron frequency v c . But simple 
models of relativistic jets in blazars suggest that just the opposite is the case (see §3.2). 
In this spirit, we propose that high-T B radiation could result from electron cyclotron maser 
instability, which operates efficiently when v c ^> u p and which has been applied successfully 
to explain auroral kilometric radiation, Jovian decametric radiation, and possibly certain 
stellar flares (Melrose 1999). 

In §2 we describe how small-scale, current-carrying magnetic mirrors, which might re- 
sult from hydromagnetic instabilities, shocks and/or turbulence, set up ideal conditions for 
electron cyclotron maser emission. As current propagates into a mirror, the motion of the 
current-carrying electrons along the flux tube is inhibited. To maintain the current through 
the mirror, a parallel electric field is established. The motion of accelerated electrons along 
converging flux tubes automatically creates the population inversion required for a cyclotron 
maser to operate. We then show that conditions thought to exist in blazar jets are favorable 
for this process to operate (§3). We estimate the fraction of the jet that needs to be filled 
with masers to produce a given observed brightness temperature, and show that even modest 
relativistic beaming (with a bulk Lorentz factor < 10) greatly increases the efficiency of the 
mechanism. Of the several processes that can inhibit the escape of high-T b maser emission, 
synchrotron absorption is potentially the most serious. In §3.2 we show that it is likely to 
prevent most of the maser radiation from reaching us, but argue that high brightness temper- 
atures can still be produced outside the synchrotron photosphere, in a thin boundary layer 
along the interface between the jet and the ambient medium. Induced Compton scattering 
provides an additional constraint on the maximum brightness temperature of a masing jet, 
while moderate relativistic beaming greatly facilitates the escape of the radiation. Finally 
(§4), we summarize our main results, generalize the model to the production of maser emis- 
sion at frequencies outside the GHz range, and discuss how the maser emission could scale 
with the mass of the system. We find that Galactic microquasars are particularly promising 
sources of cyclotron maser radiation, at frequencies that could range upward to the infrared 
or optical range. 
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2. Electron Cyclotron Maser Emission 

2.1. General conditions for cyclotron masers 

Cyclotron maser emission is thought to be produced by mildly relativistic electrons at 
frequencies close to the cyclotron fundamental or its low harmonics. Relativistic corrections 
to the gyrofrequency are essential to the operation of the maser instability (Wu & Lee 1979, 
Holman, Eichler & Kundu 1980, Melrose & Dulk 1982), which loses efficiency if the electron 
energy is too low (< 10s of keV). Masers can operate at relativistic energies (Louarn, Le 
Queau & Roux 1986), but it is possible that the effect weakens if the electron are too 
relativistic, since the resonances can become too broad. Efficient cyclotron maser emission 
also requires a highly magnetized plasma, in the sense that the electron plasma frequency, 
u p oc nj 2 , is much smaller than the cyclotron frequency, v c . Emission occurs at close to the 
cyclotron frequency if the condition v v <C v c is met. 

Several schemes have been proposed to obtain the kind of population inversion required 
for a cyclotron maser. The most common approach, called the "loss-cone" mechanism, is to 
suppose that a distribution of electrons with an initially isotropic pitch-angle distribution 
propagates down converging magnetic field lines toward an absorbing boundary, such as a 
stellar or planetary atmosphere (Holman et al. 1980, Melrose & Dulk 1982, Willes & Wu 
2004). The particles with sufficiently low pitch angles reach the boundary and are absorbed, 
while those with pitch angles above a certain threshold are reflected by the magnetic mirror. 
Loss-cone masers have the advantage that they take an initially isotropic electron distribu- 
tion function and naturally turn it into an inverted distribution through the selective loss 
of particles in a fixed magnetic field geometry. However, loss-cone inversions are readily 
quenched by pitch-angle diffusion, caused by the radiation process. Since the loss cone can 
be filled in without strong radiative losses, loss-cone masers tend to be inefficient (Melrose 
2002). 

Several authors have recently focused on the "shell" or "horseshoe" mechanism, in which 
the electron distribution is driven away from isotropy by a coherent acceleration mechanism 
(Pritchett 1984, Winglee & Pritchett 1986, Louarn et al. 1990, Delory et al. 1998, Pritchett 
et al. 1999, Bingham & Cairns 2000, Bingham et al. 2003). Since the low-energy electrons are 
depleted by the acceleration process, the distribution function resembles a shell or horseshoe 
on a v || vs. v± plot (see Fig. 1). This mechanism is more robust than the loss-cone mechanism 
(Pritchett et al. 1999), because filling in the shell (and thus quenching the inversion) cannot 
occur without strong energy losses. The maser instability growth rates tend to be higher and 
the inverted distribution can be pumped vigorously and continuously by a parallel electric 
field driving a current into a mirroring field geometry. Ergun et al. (2000) found strong 
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evidence from in situ observations that auroral kilometric radiation results from a shell 
maser, rather than a loss-cone maser. No absorbing boundary is needed for a shell maser, 
which is an advantage for our application since such a boundary is unlikely to exist in the 
low-density regions of a relativistic jet. 

2.2. Cyclotron maser emission from current-carrying magnetic mirrors 

We suppose that a blazar jet contains strongly magnetized, highly turbulent plasma. 
The turbulence could be driven by boundary conditions at the base of the jet (e.g., flux loops 
erupting from an accretion disk), hydromagnetic instabilities, strong pressure gradients, or 
velocity shear internal to the jet or near the jet boundary. The result is a highly tangled 
magnetic field topology containing strong, mostly field-aligned currents, current sheets, and 
flux ropes such as seen in the interstellar medium (e.g., Spangler 1999) and in active regions 
in the solar corona (e. g., Dmitruk, Gomez & Deluca 1998). These currents inevitably local- 
ize and dissipate in relatively confined regions (Parker 1988). As a result, large contrasts in 
magnetic field strength and direction are inescapable. A change in magnetic field strength 
along a flux rope sets up a magnetic mirror. Field-aligned currents in magnetic mirrors are 
expected in regions of oblique shocks (Bingham et al. 2003), magnetic field reversals asso- 
ciated with current sheets, strong pressure gradients, boundaries and shear layers, magnetic 
reconnection, or instabilities (kink or sausage) related to quasi-force-free flux ropes. Such 
magnetic field topologies, for example, are often observed in the solar wind. Magnetic field 
measurements in the solar wind during active periods (e.g. Tsurutani et al. 1995; Lepping et 
al. 1997) display order of magnitude variations in |B| and rapid (or small-scale) changes in 
magnetic field direction. The large variations in the magnetic fields are often associated with 
erupting flux ropes (e.g., Lepping et al. 1997) or shocks in corotating interaction regions 
(Tsurutani et al. 1995). It is a basic premise of this article that current-carrying mirrors 
are fundamental to strongly turbulent plasmas but the dominant mechanism that creates 
the field-aligned currents and changes in |B|, while interesting and important, is for future 
debate. 

Below, we demonstrate that a field-aligned current in a magnetic mirror creates ideal 
conditions for a shell maser. The current-carrying magnetic mirror sets up a parallel electric 
field and causes the accelerated electron distribution to evolve into a shell in velocity-space 
which has been shown to produce powerful maser emissions. In a strong magnetic mirror 
(mirror ratio R > 10), such as in the stellar or planetary cases, the parallel potential is large 
resulting in strong electron acceleration. Furthermore, the electron distribution can fully 
evolve into an unstable shell distribution resulting in a powerful maser. A weak magnetic 
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mirror (R < 2), however, leads to weaker electron energization and the electron distribution 
does not fully evolve into a shell, so the maser is far less powerful. In a strongly turbulent 
plasma, the current within the flux ropes or current sheets will encounter a wide range 
of magnetic mirrors, the vast majority of which will be weak. We show, however, that a 
relatively sparse packing of moderate magnetic mirrors (R ~ 5) can produce an extremely 
bright (Tb ~ 10 15 K) source. To quantify the problem, we focus on the simple case of a 
slender, quasi-force-free flux rope characterized by a transverse (toroidal) magnetic field B 
and a radius r that has a transverse (toroidal) mirroring field B m with R = B m /B ~ 5. 
A cartoon depicting the flux rope is shown in Figure 1. If the mirror has a characteristic 
radius r m , the longitudinal current density required to maintain this external field is 

jo ~ outside the mirror; and j m ~ — — inside the mirror. (1) 

4tt r 47r r m 

In a mildly relativistic (non-positron) plasma, electrons carry the majority of the field- 
aligned currents. As the flux rope evolves, the mirror force restricts the available phase space 
in the electron distribution outside of the mirror region that can carry currents through the 
region of high magnetic field (Fridman and Lemaire 1980). As in the case of Earth (e.g., 
Weimer et al. 1987, Lyons 1980, Elphic et al. 1998) and Jupiter (e.g. Su et al, 2003 and 
references therein), a parallel electric field will develop if the mirror ratio is such that 

c B n p e ec 

Jo = ~ A > — b - > y 2 ) 

An r R 

where n is the electron density in the flux rope and f3 e c xs {kTf.jm^l 2 is the electron thermal 
speed. Since strong maser activity requires mildly relativistic electrons, we will assume 
electron temperatures T e > 100 keV. Such ambient temperatures could be maintained by 
dissipation at the current sheets, given that the radiative cooling time scales are typically 
longer than flow times under blazar conditions. If R < {mi/rrieY^iTe/Ti) 1 / 2 (where subscript 
% refers to the ions), the electron drift speed exceeds the ion thermal speed. Such currents 
are ubiquitous in the Earth's auroral zone and magnetosphere (e.g., McFadden, Carlson & 
Ergun 1999). Electron drift speeds routinely are observed to exceed ion thermal speeds and 
approach, but rarely exceed, the electron thermal speed. In observed regions of electron 
cyclotron maser activity, these high drift speeds excite a number of plasma modes including 
acoustic turbulence, ion cyclotron waves, lower hybrid waves, and quasi-electrostatic whistler 
emissions (see Kindel & Kennel 1970 for a review). These emissions, however, do not offer 
sufficient resistivity to disrupt the electron current nor do they hinder the electron cyclotron 
maser process. 

The electric potential $ along the magnetic field required to maintain the current 
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through the pinch can be expressed as 



Rj = , ® e ] n \, lo with \<^-<R (3) 
J ° (2nm e kT e y/ 2 ~ kT e ~ K ' 



(Knight 1973, Fridman and Lemaire 1980) which, when combined with eq. (1), becomes 

V2 / - x 1/2 



B m e$ ( e$ \ 1 ( 2e$ 



(- 



Alienor m {2-n) 1 / 2 j3 e m e c 2 \4nkT e J \m e c 2 



(4) 



The parallel electric field increases the current density to the required level by accelerating 
electrons along the magnetic field lines. The right hand side of eq. (4) is broken into two 
factors. The second factor is simply the electron drift velocity from acceleration by the 
parallel potential. The leading factor is a correction due to the mirroring geometry which 
slows the electrons as they move into regions of high magnetic field. 

The above equations are valid in the subrelativistic limit. The growth rate of the maser 
instability increases steeply with electron energy up to relativistic energies, and levels off or 
declines at higher energies. If the initial electron energies are already mildly relativistic (our 
reason for assuming T e > 100 keV), then acceleration up to electron Lorentz factors 7 > 2 
requires modest mirror ratios, R ~ 5. The maser instability is so powerful that whenever it 
operates, electrons entering the pinch are unlikely to exceed mildly relativistic energies and 
we will use the nonrelativistic formulae throughout. We will therefore assume, for purposes of 
quantitative estimates below, that e<3> ~ 5kT e ~ m e c 2 and -B m /(47ren r m ) = j m /(ecn ) ~ 1. 

Because the increase in energy is the same for all electrons, the Maxwellian velocity 
distribution is not preserved in this acceleration process. In addition to the characteristic 
energy of the distribution shifting to higher values, the low-energy electron population is de- 
pleted, leading to a situation in which df{y)/dv\\ > 0. Although this represents a population 
inversion, it is not one that favors cyclotron maser instability. In order to obtain a cyclotron 
maser a second step is required in the evolution of the electron distribution function — con- 
version of the inversion in vn to one in v± (i.e., df(v)/dv± > 0). This occurs naturally if the 
accelerated electron current propagates into a region of sufficiently increasing magnetic field 
strength: a magnetic mirror. 

Figure 1 illustrates the population inversion process. The parallel electric field forms on 
the low-potential side of the pinch region as on Earth and Jupiter. The unaccelerated electron 
distribution has a small loss cone which is weakly unstable to the maser (Fig. lb). The 
accelerated electron distribution rapidly evolves in pitch-angle, conserving energy {m 2 c A + 
p 2 c 2 ) x / 2 and magnetic moment /i = p\/ (me|B|), where p is the momentum. Near the highest 
field region of the mirror, the distribution resembles a horseshoe in 2-D or a hollow shell with 
a loss cone in 3-D (Fig. lc). The resonance condition of the maser (dashed line in Fig. lc) 
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is for purely X-mode emission with a frequency v = v c /l- The shell maser draws its energy 
from the entire inverted electron population and can therefore convert a substantial fraction 
of the j • E power into electromagnetic waves. 

Both the acceleration and mirroring are likely to occur in close proximity or co-spatially. 
Once a strong vj_— inversion occurs, the stimulated emission of cyclotron waves runs away 
on a very short timescale and quickly saps the available free energy, filling in the low-energy 
portion of the electron distribution and essentially thermalizing the electrons as they pass 
through the mirror region. However, this does not quench the maser: the shell maser process 
is continuous (Ergun et al. 2000b) with the maser output proportional to the supply of j • E 
work associated with particle acceleration through the mirror. 



The emission from a single cyclotron maser is highly anisotropic. In the case of a shell 
maser it is strongly concentrated in the plane normal to the local magnetic field. While the 
spectrum peaks at just below the cyclotron frequency, the bandwidth Au/u ~ (7 — l)/7 
(where 7 is the random electron Lorentz factor) is as large as 1/2 if the radiating electrons 
are moderately relativistic (7 ~ 2). In practice, the detailed emission properties of any one 
maser site is unimportant, because the observed radiation will be the incoherent sum of 
emission from an extremely large number of masers. The angular distribution of radiation 
will be largely isotropized by the randomness of the magnetic field directions and the overall 
bandwidth of coherent radiation will be determined primarily by the range of magnetic field 
strengths. A randomized distribution of magnetic field directions would also largely wash out 
the strong circular polarization associated with a single cyclotron maser, although any bias 
in the mean field properties (e.g., a preferred direction relative to the line of sight) should 
show up in the polarization. In the following, therefore, we will estimate the properties of a 
typical maser assuming that it radiates isotropically. 

The power of a single maser is a fraction £ of the power dissipated in the mirror, 



where we have assumed that the acceleration volume is cylindrical with radius r m . Based 
on our discussion above, we estimate <3> ~ m e c 2 /e and adopt j m ~ cB m / (47rr m ) ~ n ec. The 
efficiency of a shell maser, £, is limited thermodynamically to £ max ~ 0.6(e$ — /cT e )/e<3>, but 
this thermodynamic limit is never reached in practice. The maser efficiency can be reduced 
by competing modes that draw their energy from the unstable shell distribution. However, 
almost all of the competing modes are at frequencies below the electron plasma frequency, so 



2.3. Maser emissivity 




(5) 
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the maser condition v v <C v c (Melrose & Dulk 1982) ensures that growth rates of competing 
modes are small. In the case of the Earth's aurora, with v v jv c ~ 10~ 2 , satellite measurements 
indicate that the maser efficiency can be as high as a ~ 1 per cent (Gurnett 1974). In the 
case of blazar jets, v v j v c could be even smaller (see §3.2) and the efficiency could be as large 
as > 10%. We will adopt the normalization £_ 2 = £/0.01. 

Expressing n in units of cm -3 and setting B 2 = B m /10 2 G, we obtain a pinch radius of 

r m ~ 1.7 x 10 10 S 2 no 1 cm (6) 

and a power per maser of 

P rn ~ ir£n m e c 3 r 2 m ~ 2.1 x 10 23 £_ 2 5 2 n 1 erg s" 1 . (7) 

Although our immediate aim is to explain masers observed at frequencies of a few GHz, we 
chose a fiducial pinch magnetic field strength of 100 G, corresponding to a cyclotron frequency 
~ 0.3 GHz, because the maser emission from a blazar jet is expected to be Doppler boosted 
by a factor ~ 10. If the maser is modeled as a sphere of radius r m , the mean intensity is 

h ~ P 2 m 2 ~ 1.7. x 10- 8 ^_2 J B 2 - 1 n erg cm" 2 s" 1 Hz^ 1 sr" 1 . (8) 

107T V c T m 

The associated brightness temperature is 

T B = ~ 6.9 x 10 11 e_ 2J B 2 - 3 n K. (9) 

The observed intrinsic brightness temperature is an average over the unresolved source, given 
by 

(T B )obs = C m TT B , (10) 

where T is the bulk Lorentz factor of the masers toward the observer (more precisely, the 
Doppler factor) and C m is the "covering factor" of masers. Physically, this is the fraction of 
the solid angle of the source subtended by masers if C m < 1 and the mean number of masers 
along a line of sight if C m > 1. 



3. Application to blazar jets 

3.1. Relation to jet properties 

Observations of diffractive scintillation of the blazar J1819+3845 (Macquart & de Bruyn 
2004) appear to require intrinsic brightness temperatures (T B ) ohs > 10 14 K at frequencies of a 
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few GHz. This estimate comes from an upper limit on the size of the scintillating component 
< 10 uas; the actual brightness temperature could be far higher. 

We assume that this emission comes from cyclotron masers in a jet moving toward us 
with a bulk Lorentz factor T = 10IY The observed frequency of the maser emission is 
therefore z/ obs = t'GHz GHz ~ Yv c ~ 2.8TiB 2 GHz. Because the maser emission is produced 
in magnetic mirrors with R ~ 5, the maser frequency can be much larger than the (Doppler- 
boosted) cyclotron frequency associated with the mean magnetic field strength B in the flux 
rope. For the purpose of our illustrative numerical examples below, we identify the latter 
with the mean magnetic field in the jet. However, one should bear in mind that some maser 
emission could come from mirrors associated with current sheets where B is lower than the 
mean. The density in the mirror, which we have taken to be n , is probably of the same 
order as the mean jet density. Using the relations given in §2.3, we can relate the properties 
of maser emission to the mean conditions in the jet. 

It is instructive to express the local conditions, n (r) and B (r), in terms of the energetics 
of the blazar jet. We assume that the jet subtends a solid angle f2 = 0.1f2_i sr as it propagates 
away from the black hole and parameterize the distance from the black hole r in terms of 
the gravitational radius of the black hole, r = xGM/c 2 = 1.5 x 10 13 M 8 x cm, where M 8 is 
the black hole mass in units of 1O 8 M . The total jet power Lj can be divided into two parts. 
The kinetic energy flux is L^e — i]nQmiC 3 T 2 r 2 fl = Lj/ (1 + a). Here, a is the magnetization 
parameter and r\ equals 1/rrii times the mean mass per electron [= im e jm/) + (nj/no), where 
rii is the ion density]. The parameter rj ranges from 1 for normal ion-electron plasma to 
m e /mi if pairs outnumber ions by more than 2000-to-one. The electromagnetic (Poynting) 
flux is L EM = (Bl/A-n)cT 2 r 2 Vt = oLj/{\ + a). Letting L 45 = Z^/IO 45 erg s -1 , we obtain 

n = 9.9 x 10 9 (1 + ay^L^T^M^nzlx- 2 cm" 3 (11) 

and 

B = 1.1 x 10V/ 2 (1 + ay^L^T^M^nZ^x- 1 G. (12) 

One may readily check whether the plasma frequency is much smaller than the electron 
cyclotron frequency. Under mean jet conditions, 

{^j = 3.2 x lO" 3 ^ ~ 2.8 x lO-V^V 172 - (13) 

Note that this ratio does not depend explicitly on r. The ratio is even smaller, by a factor 
~ R~ l , in masing regions where B m ~ RB . Thus, the maser condition is satisfied at all radii 
provided that arj 3> 2.4 x 10 3 i? 2 . (Note also that ion cyclotron damping need not quench 
the turbulence on scales larger than r , since the ion gyroradius is likely to be smaller than 
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r for the parameters under consideration, r,i/r ~ {m e T i lmiT e ) l l' 1 <j 1 rj 1 [E. G. Zweibel, 
private communication].) 

The observed frequency of the maser emission is related to the distance of the masers 
from the black hole and the characteristic mirror ratio R through 

B ~ ^^R'^ 1 G. (14) 

Eliminating B between equations (12) and (14), we find that maser emission at a frequency 
v is produced at 

r obs ~ 4.7 x lO'V^l + ay^Rv^Lli 2 ^ 2 cm, (15) 

corresponding to x obs ~ 10 3 gravitational radii for R ~ 5 and with M§ and all other param- 
eters set equal to one. 

The typical jet density at r obs is 

n (r obs ) = 9.4 x lO^rj^u^R- 2 ^ 2 cm" 3 (16) 
and we can calculate the size, brightness temperature and power of a typical maser: 

r m (r obs ) ~ 6.5 x l^ar ] v^R 2 V 1 cm (17) 

T B (r ohs ) ~ 1.5 x 10 18 £_ 2( 7-V^GH z ^~ 2r i K ( 18 ) 
P m (r ob s) ~ 2.9 x 10 17 ^ 2 a V R 2 erg s" 1 . (19) 

In terms of the observationally inferred brightness temperature T 15 = (T B ) o bs/^0 15 K, 
we find from eq. (10) that the required covering factor of masers is 

C m (r obs ) ~ 6.7 x 10- 5 C 2 W G H z i? 2 T 15 i7 2 . (20) 

The masers need cover only a small fraction of the source on the sky, because each one has 
such a high brightness temperature. Another useful quantity is the volume filling factor. 
Assuming that the depth of each maser along the line of sight is ~ r m , the volume filling 
factor of masers is f m ~ (r m /r)C m , where r is the depth of the region containing the masers. 
We presume this depth to be comparable to the distance of the masing region from the black 
hole at the base of the jet. Setting r ~ r obs , we have 

/ m (r obs ) ~ 9.2 x 10- 16 C 2 V 3/2 (1 + ( 7) 1 / 2 r ? 2 z/ GHzJ R 3 T 15 L 4 - 1/2 rr 1 fi V 1 2 . (21) 

The total number of maser sites required is 

N m (r obs ) ~ fm (y^j (^j ~ 3.0 x 10 15 C2(1 + arV^GUzT^L^ 4 . (22) 
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Therefore, the total, broadband maser power emitted in the comoving frame of the jet (at 

r ~ r obs ) is 

^em(robs) ~ N m (r ohs )P m (r ohs ) ~ 8.7 x 10 3 V(1 + ay^^R^L^ 4 erg s" 1 , (23) 

as one can verify directly from (T B ) ohs and r obs . The broadband "isotropic power" inferred 
by an observer is boosted by a factor T 4 . These energy requirements are modest for the 
conditions thought to exist in blazar jets: for the model described here the volume-averaged 
emissivity of maser radiation at r obs is only a fraction ~ 2.5 x 10~ 10 T 15 i? 2 IY 2 of the comoving 
magnetic energy density passing through r obs per dynamical time. 



3.2. Can cyclotron maser radiation escape from blazar jets? 

High-Tg radiation at ~ GHz frequencies, and cyclotron maser emission in particular, 
is subject to at least four processes that can inhibit the escape of radiation: 1) induced 
(stimulated) Compton scattering, 2) stimulated Raman scattering, 3) cyclotron absorption in 
the second and higher harmonics, and 4) synchrotron absorption by ultrarelativistic electrons. 
Synchrotron absorption is by far the most serious impediment, and will lead us to consider 
emission from a thin boundary layer. 



3.2.1. Induced Compton and stimulated Raman scattering 

Once cyclotron maser radiation is produced it must traverse the gas in the non-masing 
regions along the line of sight, in order to escape to the observer. Let us first consider the 
constraints imposed by induced Compton scattering and stimulated Raman scattering. We 
assume that we are looking down the jet, and that the main column density of intervening 
material is located close to the maser emission region, at r obs . Induced Compton scattering 
will limit the mean brightness temperature to a value that satisfies the condition 

"'{-1 B/obs -p_2 <f -1 /<->/l\ 

n ra T — r < 1, (24) 

m e c z 

where o~t is the Thomson cross-section (Wilson 1982; de Kool & Begelman 1989; Coppi, 
Blandford & Rees 1993; Sincell & Krolik 1994). In this expression, one power of T comes 
from the transformation of the observed brightness temperature to the comoving frame, and 
the other comes from the Lorentz contraction of the radial scale length r (or, equivalently, 
from the Doppler-shifted scattering rate). Note that the relevant brightness temperature 
to use in evaluating the stimulated scattering rate is the angle-averaged mean brightness 
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temperature in the comoving jet frame, rather than the higher T B associated with each 
maser beam. This is because the scattering rate is proportional to the occupation number 
of available final states for the scattered photon, which includes an angular filling factor. 
Substituting n from above we obtain an upper limit to the observed brightness temperature 
escaping from a given radius: 

(T B ) ohs < 6.0 x 10 12 (1 + a^L^TiMstt^x K. (25) 

Note the importance of relativistic beaming effects in facilitating the escape of radiation: a 
factor of 10 increase in T increases the maximum brightness temperature by a factor of 10 4 . 
If we set x = x bs from eq. (15), we obtain 

(T B ) ohs < 1.9 x lO'V^l + af/^^RL-^Vl^ll K. (26) 

Thus, induced Compton scattering can pose a serious threat to the escape of cyclotron maser 
radiation in our picture, if (T B ) ^, S is much higher than ~ 10 15 K. However, the difficulties are 
not insurmountable. As noted above, the effects of induced Compton scattering are greatly 
reduced if V is even slightly higher; escape of radiation is also facilitated by larger R, larger 
ratio of Poynting flux to kinetic energy flux o (for fixed jet power), and larger jet solid angle 
Q. Pair-dominated jets (i.e., with 77 ^> 1) have more trouble producing high (T B ) ohs , because 
of the larger electron and positron column densities required to carry a given Lke- 

Induced Compton scattering can also heat the electrons to temperatures well above the 
inverse Compton temperature associated with the total ambient radiation field (Levich & 
Sunyaev 1970). However, for the conditions we are considering one can show that ordinary 
cyclotron cooling exceeds induced Compton heating for all temperatures of interest. More- 
over, cyclotron cooling itself is not very important near r obs , since the cooling timescale 
exceeds the flow time. Our model assumes that the mean temperature of electrons entering 
the pinch is mildly relativistic (to obtain reasonable maser efficiencies); such temperatures 
would have to be maintained by dissipative processes in the flow, such as reconnection or 
shocks. 

We next consider stimulated Raman scattering, which was studied in detail by Levinson 
& Blandford (1995). As shown by their Fig. 6, Raman scattering becomes important relative 
to induced Compton scattering at low electron temperatures and high densities. (We note 
that Levinson & Blandford considered Langmuir turbulence in a non-magnetized plasma, 
whereas we have the opposite situation, u c ^> u p . Presumably, this difference between our 
case and theirs would tend to diminish the importance of Raman scattering.) At x Q b s , the 
mean electron density is n (x o bs) ~ 10 5 °"~ 1 ^~ lz/ GHz-^ > ~ 2 rr 2 cm -3 . As shown by Blandford & 
Levinson's Fig. 6, at such densities stimulated Raman scattering is completely unimportant 
compared to induced Compton scattering, for all values of {T B ) ohs and T e under consideration. 
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3.2.2. Local cyclotron absorption 

A major difficulty in applying the cyclotron maser to astrophysical radio sources is that 
emissions at the fundamental frequency are expected to be nearly entirely absorbed as they 
propagate through regions where they are at higher harmonics of the local electron gyrofre- 
quency (Melrose & Dulk 1982). There are several possible windows of escape: nearly parallel 
propagation (R-mode) where the optical depth is small, nearly perpendicular propagation 
with partial mode conversion, or mode conversion and other nonlinear processes in or near 
the absorption layer (Robinson 1989). We suggest that cyclotron absorption may be avoid- 
able in current-carrying mirrors because the thickness of the second harmonic absorption 
layer is small. A shell maser produces emission in the X-mode propagating perpendicular 
to B. As the radiation moves away from the masing region, the magnetic field falls roughly 
as B = B m r m /r. Correspondingly, the maser emission travels through the second harmonic 
layer at r = 2r m . The damping rate {vd) at this layer is unlikely reach as high as u p , which 
in turn is much smaller than v c . The thickness of the second harmonic damping layer can 
be estimated as ~ r m u p /u c so the travel time through the region is r = r m u p /u c c yielding an 
attenuation coefficient of ~ r m v^jv c c. 

From equation (6), the size of the mirror can be expressed as r m ~ O.lcz^c/Vp, so the 
attenuation coefficient is roughly 0.1. Thus, magnetic mirrors along current-carrying flux 
ropes not only form natural environments for the shell maser inversion process but also allow 
for the escape of the radiation because the regions of enhanced field are small and the field 
drops rapidly with transverse distance, resulting in a thin second harmonic absorbing layer. 

3.2.3. Synchrotron absorption 

If the magnetic field in the masing region is only a few times higher than the mean field 
outside, as we have argued, then there is no clearcut distinction between cyclotron absorption 
in the first few harmonics (where the absorbing electrons have 7 > 1) and synchrotron 
absorption (where 7 ^> 1) at radii r > r G b s - At radii ^> r G b s along the line of sight, where the 
magnetic field has dropped by a large factor, the dominant absorption will be squarely in 
the synchrotron regime, and we will make a relatively small error by using the synchrotron 
formulae everywhere. 

The flat radio spectra exhibited by many blazar jets are generally assumed to arise from 
self-absorbed synchrotron radiation. Typical intrinsic brightness temperatures < 10 11 K at a 
few GHz (Readhead 1994) indicate that the emission is produced where the mean magnetic 
field -Bq is considerably lower than in the masing regions, and therefore at distances r ^> r Q b s 
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from the black hole. If this emission is truly self-absorbed and the emitting regions lie along 
the line of sight to the masers, then they will block the escape of the maser emission. 

Therefore, it should come as no surprise that GHz maser radiation can only escape from 
the interior of a jet with an anomalously low synchrotron emissivity (compared to the usual 
assumptions), corresponding to a very low ratio of relativistic electron energy density U rc \ 
to magnetic energy density Ub = Bq/8tt. Recall that standard models for self-absorbed jets 
(Readhead 1994) assume that U TC \ and Ub are close to equipartition. If U rc \ <C Ub, then 
the synchrotron emissivity saturates at a lower brightness temperature and the synchrotron 
photosphere lies closer to the black hole along the jet axis. (These arguments apply to the 
intrinsic brightness temperature; the perceived brightness temperature also depends on the 
Doppler factor, which could vary from jet to jet.) If the maser emission from the jet interior 
is blocked by synchrotron absorption, then we will see only the small fraction of emission 
produced outside the jet photosphere, presumably from a thin boundary layer along the wall 
of the jet at radii ~ r obs . 

To estimate the energy density of relativistic electrons consistent with the escape of 
maser emission from the jet interior, let us suppose that the relativistic electrons have a 
power-law distribution in random Lorentz factor n re i(7) oc 7~ 2 for 7 min < 7 < 7max- If 7min 
is small enough that the corresponding synchrotron critical frequency is below the observing 
frequency (as Doppler-shifted into the comoving frame), 7 m i n < 7 C (^, -Bo) — 5-B 1 ^ 2 i / gbz^i ^ 
(corresponding to 7 min < r/r b s ), then the synchrotron absorption coefficient at the observing 
frequency z/qhz GHz is given by (Rybicki & Lightman 1979) 

-^^^R^t^--" 1 (27) 

In the magnetically dominated jet environment we might reasonably expect the energy den- 
sity of relativistic electrons not to exceed the thermal energy density, Uth ~ 3riokT e /2. 
Normalizing U TC \ to U t h instead of Ub we obtain the optical depth along the line of sight to 
radius r > r obs , 

(\ 7 8.1 n9 ^ 1/2 (l + tr)- 1/2 -1 U Tcl kT e 3 , 1/2 1/2 fry 3 
Ts y nch(r) ~ 7.8 x 10 —— —r] — — — — R r\ i2_/ . (28) 

hl(7max/7min) U th 100 keV Vjobs/ 

If 7min > 7o both the absorption coefficient and the optical depth are lowered by a 
factor (7min/7c) 8 ^ 3 ) implying that even a moderate increase in 7 min can greatly suppress the 
absorption of maser radiation. Nevertheless, it is clear that synchrotron absorption provides 
by far the most serious impediment to the escape of maser radiation. To have any hope of 
observing high brightness temperature emission, we must consider rather "optimistic" values 
of the main parameters. One might reasonably assume R ~ 5 (although higher values are 
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not excluded), ln(7 max /7 min ) ~ 10, and r] ~ 1 (electron- ion plasma). We shall take kT e ~ 100 
keV and — more speculatively — suppose that U re \ could be as low as < 10~ 3 U t h, which is 
consistent with a Maxwellian tail on the thermal electron distribution (although Compton 
cooling on diffuse ambient radiation could deplete such a tail). The optical depth is reduced 
to 

6.2x10 U IC \ „_ lr l/2 n -l/2/ 1 . 8/3 \-l / on \ 

r sy nch(r bs) ~ a i/ 2 ( 1 + (T y/2 W*U^ 1 45 ( 7min) ' ( } 

Note that r sync h(r obs ) is independent of frequency, since the decrease of column density with 
radius is offset by the decreasing frequency of the emission. 

There are combinations of parameters — large values of a, T and 7 m i n ; or smaller values 
of U re \ — that lead to r sync h < 1, thus allowing the maser emission to escape from the interior 
of the jet. None of these quantities (with the possible exception of T) is even moderately 
well-constrained, especially at radii ~ r G b s that are smaller than the synchrotron emission 
radius. Moreover, bright maser emission may still be observable even if r syIlc h > 1. For 
example, the maser emission could escape from a thin boundary layer along the wall of the 
jet. The fraction of the jet volume involved in producing observable maser radiation at r 
would then be ~ 1/T sync h(r); the filling factor of maser sites in the boundary layer would have 
to be correspondingly higher. Indeed, preferential formation of maser-producing mirrors near 
the jet boundary is not implausible, since strong shear, tending to stretch and compress the 
field, is likely to be especially prevalent near the walls of the jet. An attractive feature of the 
boundary layer model is that shear at the interface with the ambient medium could drive 
supermagnetosonic turbulence, leading to localized magnetic field strengths far higher than 
the mean jet field. This would allow the observed maser radiation to be produced at larger 
radii, facilitating its escape. 

While it would certainly help for maser emissivity to be enhanced close to the walls of 
the jet, compared to the interior, it is by no means necessary. This is because the maser 
mechanism is so efficient: for maser emissivity spread randomly through the jet, the ratio of 
volume-averaged maser emissivity to magnetic energy supply need only be ~ 10~ 8 Ti5lY 2 to 
provide the observed brightness temperature. If we only see the masers in a boundary layer, 
the required emissivity is increased by a factor r synch (r obs ) but is still quite small. We note 
that these estimates do not take into account the possibility of synchrotron absorption and 
induced Compton scattering in the ambient medium that surrounds the jet, the nature of 
which is very uncertain. 
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4. Discussion 

We have outlined a set of physical processes by which high-brightness temperature elec- 
tron cyclotron maser emission can be produced in blazar jets, as well as other tenuous, highly 
magnetized astrophysical media. The population inversion leading to maser emission occurs 
in transient, small-scale, current-carrying magnetic mirrors with mirror ratios > 3. The 
latter can arise from hydromagnetic instabilities, shocks, and/or turbulence driven, for ex- 
ample, by shear near the interface between the jet and the ambient medium. The mechanism 
we have described is an example of a "shell maser" , in which the population inversion results 
mainly from electron acceleration in a parallel electric field; the more commonly discussed 
"loss cone" effects are secondary. As a precondition for efficient maser emissivity, a jet should 
have a ratio of Poynting flux to kinetic energy flux that is not too small, so that the elec- 
tron cyclotron frequency exceeds the electron plasma frequency by a large factor. Another 
condition for high maser efficiency is that the thermal electrons carrying the current into 
the magnetic mirrors be hot, with kT e ~ 100 KeV. The mechanism works much better in a 
jet composed of electron-ion plasma than in a pair-dominated jet, since the latter will have 
a higher ratio of plasma frequency to cyclotron frequency and a larger induced Compton 
scattering depth. 

We are unable to estimate the volume fraction of a turbulent jet that will have prop- 
erties conducive to maser emission, but we show that GHz brightness temperatures ~ 10 15 
K can arise from systems with an extremely small filling factor of such sites. Moreover, the 
fraction of internal (magnetic) energy that needs to be converted to maser radiation per flow 
time is extremely small. These conclusions are encouraging given the severe constraints on 
the escape of maser radiation from the jet, once it is produced. Induced Compton scatter- 
ing places upper limits on the brightness temperature of the escaping radiation; however, 
the main limitation comes from synchrotron absorption by relativistic electrons within the 
jet. For all but the most extreme parameters the jet will be very optically thick to maser 
radiation, but we argue that a thin boundary layer along the wall of the jet could produce 
high brightness temperatures without excessive demands on the energetics or filling factor 
of masing regions. We note that even the boundary layer model requires relativistic electron 
energy densities far below equipartition — this could be due to weak particle acceleration 
in these inner jet regions, where little of the magnetic energy has been dissipated, possibly 
combined with the effects of rapid synchrotron or inverse Compton cooling. 

Our model is parameterized in terms of maser emission at a few GHz, but the process 
could clearly operate at other frequencies as well. According to our jet model, the radius at 
which maser radiation of a given frequency originates [r b s (V)] scales as u^ 1 . Thus, maser 
emission at frequencies of up to ~ 100 GHz could be produced near the bases of blazar 
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jets. The energy efficiency of maser emission (denned at the end of §3.1) is independent 
of u, as is the characteristic synchrotron optical depth. However, the maximum brightness 
temperature permitted by induced Compton scattering scales as u^ 1 . 

The relations derived in §3 can be scaled to relativistic jets from stellar-mass systems, 
such as microquasars. If we assume that the jet power, L, scales linearly with the mass 
of the black hole, then the frequency of maser radiation produced at a given dimensionless 
radius, x ^ s (in units of the gravitational radius GM/c 2 ), scales as v oc M -1 / 2 . Therefore, 
microquasars (and other types of X-ray binary) could emit cyclotron maser emission in the 
infrared or optical band, from the base of the jet or a corona above the accretion disk. The 
characteristic synchrotron optical depth scales as M 1//2 , independent of v\ therefore, typical 
values of r synch should be 3-4 orders of magnitude lower than those estimated for blazars. 
This greatly relaxes the constraints on boundary layer emission, and raises the possibility that 
microquasar jets may be optically thin to maser radiation from their interiors. The induced 
Compton limit on the brightness temperature emitted at a given x bs is independent of mass, 
but at a fixed observing frequency the brightness temperature limit scales as M -1 / 2 . Thus, 
high brightness temperature maser emission has an easier time escaping from lower mass 
systems. 

The principal diagnostics of planetary and stellar cyclotron maser emission — high cir- 
cular polarization and extreme variability — are probably not useful for analyzing cyclotron 
maser emission from jets. The enormous number of independent maser sites, spread over 
a region vastly larger than the size of each maser, would presumably wash out such signa- 
tures. Clues to a coherent emission mechanism could be obtained by brightness temperature 
measurements as a function of frequency as well as evidence for an unusual source geometry, 
such as a narrow boundary layer. Although the production and escape of cyclotron maser 
radiation from relativistic jets is by no means inevitable, the search for coherent emission 
remains exceedingly worthwhile because of the constraints it would impose on the structure 
and composition of these energetic flows. 

We are indebted to Ellen Zweibel, Roger Blandford and Paolo Coppi for their insightful 
comments and criticisms, and to Jean-Pierre Macquart for bringing the latest scintillation 
measurements to our attention. MCB acknowledges support from National Science Founda- 
tion grant AST-0307502. 
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Fig. 1. — Fig. 1. (a) Cartoon of a current-carrying magnetic mirror on a quasi-force-free 
flux rope. The nearly field-aligned current increases both the toroidal and axial components 
of the magnetic field. Electrons carry the current of the flux rope from top to bottom. A 
parallel electric field, as in the aurorae on Earth and Jupiter, is required to maintain the 
electron flow through the mirror region. The combination of parallel potential and magnetic 
mirror evolves the initial electron distribution (b) into a horseshoe-shaped distribution (c) 
when viewed in 2-D (a shell when viewed in 3-D). In these plots we use a mirror ratio R = 5, 
a current \ J zm \ = 30mA/ m 2 (\J z o\ — 6m A /m 2 ), and a 500 keV parallel potential which is 
consistent with the current for an electron temperature of ~ 100 keV and density of 100 
cm" 3 . The adiabatic, static Vlasov code used for these calculations is described by Ergun 
et al. (2000b). 



